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Fast modulation and switching of light at visible and near-infrared (vis-NIR) frequen¬ 
cies is of utmost importance for optical signal processing and sensing technologies. No 
fundamental limit appears to prevent us from designing wavelength-sized devices ca¬ 
pable of controlling the light phase and intensity at gigaherts (and even terahertz) 
speeds in those spectral ranges. However, this problem remains largely unsolved, 
despite recent advances in the use of quantum wells and phase-change materials for 
that purpose. Here, we explore an alternative solution based upon the remarkable 
electro-optical properties of graphene. In particular, we predict unity-order changes 
in the transmission and absorption of vis-NIR light produced upon electrical doping of 
graphene sheets coupled to realistically engineered optical cavities. The light intensity 
is enhanced at the graphene plane, and so is its absorption, which can be switched 
and modulated via Pauli blocking through varying the level of doping. Specifically, 
we explore dielectric planar cavities operating under either tunneling or Fabry-Perot 
resonant transmission conditions, as well as Mie modes in silicon nanospheres and lat¬ 
tice resonances in metal particle arrays. Our simulations reveal absolute variations in 
transmission exceeding 90% as well as an extinction ratio > 15 dB with small insertion 
losses using feasible material parameters, thus supporting the application of graphene 
in fast electro-optics at vis-NIR frequencies. 


I. INTRODUCTION 

Graphene -the two-dimensional (2D) honeycomb lat¬ 
tice of carbon atoms- exhibits extraordinary optoelec¬ 
tronic properties derived from its peculiar band struc¬ 
ture of massless charge carriers [T]. Notably, its optical 
absorption can be switched on/off via electrical doping. 
In its undoped state it absorbs a fraction ira « 2.3% 
of the incident light Hug over a broad spectral range 
within the vis-NIR as a result of direct electron-hole pair 
transitions between its lower occupied Dirac cones and 
the upper unoccupied cones (two inequivalent ones in ev¬ 
ery Brillouin zone [US]). In contrast, when electrically 
doped, an optical gap is opened that suppresses verti¬ 
cal optical transitions for photon energies below 2|FJ?|, 
where E-p is the change in Fermi energy relative to the 
undoped state (see Fig. Op. In practice, values of Ep 
as high as 1 eV can be obtained through electrical gating 
[5], therefore enabling the modulation of light absorption 
down to the visible regime. Chemical methods permit 
achieving even higher levels of doping [6., which could be 
combined with additional electrostatically induced vari¬ 
ations of Ey around a high bias point to reach control 
over shorter light wavelengths. 

Fast light modulation at vis-NIR frequencies can find 
application in optical signal processing and interconnect 
switching, where there is a great demand for integrated 
wavelength-sized devices capable of operating at tera¬ 
hertz commutation rates. The extraordinary electro- 
optical response of graphene provides a key ingredient 
for the realization of these types of devices. However, 
the exploitation of atomically thin carbon films for light 
modulation faces the problem of their relatively weak in¬ 


teraction with light. A possible solution to enhance this 
interaction is to use the intrinsic plasmons that show up 
in the optical gap of this material when it is highly doped 
urn- Resonant coupling to graphene plasmons can even 
result in complete optical absorption m, as exempli¬ 
fied by the observation of large tunable light modula¬ 
tion at mid-IR frequencies in periodically nanostructured 
graphene mm- The extension of this strategy down to 
the vis-NIR spectral domain remains a challenge, as it re¬ 
quires to laterally pattern the carbon film with <10 nm 
features, which are currently unattainable through con¬ 
ventional lithographies, although chemical self-assembly 
might offer a viable way of producing the required struc¬ 
tures [21]. 

An alternative solution consists in amplifying the ab¬ 
sorption of undoped graphene either by increasing the re¬ 
gion over which light interacts with it or by coupling the 
carbon film to an optical cavity of high quality factor (i.e., 
by trapping light during long times near the graphene). 
A broadband modulator has been demonstrated with the 
former approach by exposing a long path of an optical 
waveguide to electrically gated graphene m • Addition¬ 
ally, coupling to photonic cavities has been explored us¬ 
ing plasmonic structures, photonic-crystals, and meta¬ 
materials [231429] . For example, monolayer graphene in¬ 
tegrated with metallic metasurfaces has been used to 
control the optical response (resonance position, depth, 
and linewidth) at mid-IR frequencies [2SH23 EHj - Sim¬ 
ilarly, large intensity modulations (> 30%) of mid-IR 
light over a 600 nm bandwidth have been reported in 
graphene-loaded plasmonic antennas m- Additionally, 
a resonance wavelength shift ~ 2nm accompanied by a 
4-fold variation in reflectivity has been observed in the 
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NIR by coupling graphene to a photonic crystal cavity 
El [30]. Enhanced visible light absorption in graphene 
has also been demonstrated (without modulation) by 
combining monolayer graphene with metamaterials [31], 
gold nanovoid arrays m, and photonic waveguides m , 
as well as by coupling multilayer graphene under total 
internal reflection [33]. 

In this work, we study four different mechanisms that 
produce resonant enhancement in the absorption of un¬ 
doped graphene over the vis-NIR spectral domain, thus 
serving as optical modulators with large depth in that 
frequency range. Specifically, we focus on the coupling 
of graphene to (1) resonant tunneling transmission cav¬ 
ities, (2) resonant Fabry-Perot cavities, (3) Mie modes 
of silicon spheres, and (4) lattice resonances in periodic 
particle arrays, which we investigate by calculating reflec¬ 
tion, transmission, and absorption spectra of structures 
containing either doped or undoped graphene films. We 
predict modulation depths in vis-NIR light transmission 
exceeding 90%, with small insertion losses, thus reveal¬ 
ing the potential of graphene for fast electro-optics within 
such a technologically important range of optical frequen¬ 
cies. 


II. RESULTS AND DISCUSSION 

A. Graphene optical switch based upon resonant 
tunneling transmission 

We illustrate the concept of resonant switching and 
modulation of graphene absorption by coupling to a high- 
quality-factor planar cavity. In particular, we consider 
the multilayer structure depicted in Fig. Eh consisting of 
a high-refractive-index BN planar waveguide (tibn = 2.1) 
flanked by low-index silica spacers (nsio 2 = 1-457). The 
waveguide hosts guided modes that can be resonantly 
coupled to light of well-defined parallel wave vector (i.e., 
for a collimated incident beam). In our case, light is 
incident from the left under total internal reflection con¬ 
ditions at the BF11-Si02 interface (riBFii = 1.61). The 
evanescent spill out of light intensity penetrating inside 
the left silica spacer can reach the BN waveguide, where 
it is amplified to further extend towards the rightmost 
interface. In the absence of absorption, full transmission 
can always be achieved at a resonant wavelength that de¬ 
pends on incidence angle. This phenomenon, known as 
resonant tunneling transmission, was previously explored 
with electron waves [35] • There is complete analogy be¬ 
tween TE light propagation in the planar structure under 
consideration and the evolution of an electron according 
to Schrodinger equation [36] (see Methods). The equiva¬ 
lent electron has energy E and evolves along a potential 
profile as shown in Fig. [ 1 J 3 . The latter is directly related 
to the refractive index, with higher index correspond¬ 
ing to lower values of the potential. The presence of a 
bound state is always guaranteed in a ID cavity, and so 
is the existence of a full transmission resonance when this 


bound state lies inside the potential barrier E3- Under 
complete-transmission conditions, the intensity has to de¬ 
cay exponentially from the waveguide to the far medium 
(i.e., across the rightmost silica barrier), to reach the 
same value as the incident intensity, so that the near 
field has to be strongly amplified at the central waveg¬ 
uide. This type of enhancement, which is clearly illus¬ 
trated in Fig. [ljl, has been experimentally corroborated 
by measuring a > 100-fold increase in the fluorescence 
from quantum dots placed near the central waveguide 
under resonance conditions in a structure similar to the 
one considered here [36]. This effect can find application 
to sensing, for example by replacing the waveguide by a 
high-index fluid with dispersed analytes in it. Instead, 
in this article we use the resonant tunneling transmission 
concept to amplify the effect of absorption taking place 
at the graphene. 

The structure under consideration (Fig. Eb contains a 
graphene film on either side of the central BN waveg¬ 
uide. Besides its high index of refraction, the choice 
of BN for the central waveguide is convenient because 
this combination of materials is compatible with high- 
quality graphene m, which can be realistically described 
with the models for the conductivity a discussed in the 
Methods section. Nevertheless, we assume a conserva¬ 
tive value of the graphene mobility throughout this work, 
fi = 2000cm 2 /(Vs). When the carbon layer is highly 
doped (E-p = 1.1 eV), it becomes nearly lossless (i.e., 
small Re{cr}) at the waveguide resonance wavelength, so 
that the peak transmission reaches ~ 95% (Fig. El) and 
the light intensity enhancement at the waveguide exceeds 
a factor of 140. In contrast, in the undoped state, the car¬ 
bon layer becomes lossy (i.e., nearly real a « e 2 /Ah), so 
the enhancement is strongly suppressed, and the trans¬ 
mission drops to very small values. The extinction ratio 
(i.e., the ratio of transmissions in doped and undoped 
states) is > 15 dB. The transmission can be in fact tuned 
continuously between these two extreme values by vary¬ 
ing the level of doping (see Fig. [l^). The decrease in 
transmission produced when moving from highly doped 
to undoped graphene is due to both absorption and re¬ 
flection, as the local change in the response of the carbon 
layer produces a departure from the conditions of reso¬ 
nant tunneling. Actually, reflection accounts for the bulk 
of the depletion in transmission, as shown in the supple¬ 
mentary Fig. [6] This can be exploited to simplify the 
structure, which still undergoes unity-order modulation 
of the reflection upon graphene doping after removing 
the rightmost BF11 out-coupling medium (see supple¬ 
mentary Fig. [8]). Even more, only a single graphene layer 
is needed to modulate the structure (see supplementary 

Fig.0. 

The wavelength of operation of this modulator is es¬ 
sentially determined by the waveguide mode, as coupling 
to the BF11 media is just producing a slight shift. Under¬ 
standably, the reflection minimum is observed to be only 
mildly modified when the rightmost glass is removed (cf. 
supplementary Figs. [7] and [8|. Then, we find it useful 
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FIG. 1: Graphene optical switch based on resonant tunneling transmission, (a) Doping-induced absorption switching 
effect used in this work: we compare undoped graphene (upper scheme, Fermi level at the Dirac point), which can absorb photons 
(vertical arrow) over a broad spectral range via interband electron transitions, and doped graphene (lower scheme), in which 
Pauli exclusion blocks photon absorption when the Fermi energy Ep exceeds half the photon energy, (b) Planar multilayer 
structure considered for resonant tunneling transmission of light, including a central BN planar waveguide (not to scale) and 
two single-layer graphene films intercalated at the BN/SiCb interfaces, (c) Potential in the equivalent Schrodinger model (see 
main text and Methods), (d) Electric field intensity normalized to the external light intensity for an incidence angle of 71° 
and a free-space wavelength of 689 nm. Light is s (TE) polarized and incident from the left. Results for different levels of 
doping are offered (see legend in (e)). (e) Transmission spectra of the multilayer structure at 71° incidence for different levels of 
doping. The transmission maxima are in excellent agreement with the analytical expression offered in the Methods section (see 
arrows). The numerical labels correspond to the ratio Re{a}/(e 2 /4h) evaluated at a wavelength of 689nm. (f) Transmission 
as a function of incidence angle and wavelength for doped and undoped graphene. 


to derive a simple analytical expression for the variation 
of the waveguide resonance wavelength (see Eq. © in 
the Methods section), in which the graphene conductiv¬ 
ity enters to first order as oc Im{cr}. The role of the real 
and imaginary parts of a is thus clearly established: the 
former determines the depth of the modulation, whereas 
Im{cr} is responsible for the resonance shift. The res¬ 
onance wavelengths given by Eq. Q are indicated by 


downwards arrows in Fig. [l^, in excellent agreement with 
the observed transmission maxima. Notice the initial 
redshift with increasing doping, followed by a blueshift 
back to the original position in the perfectly transmit¬ 
ting structure, which essentially mimics the evolution of 
Im{cr} with doping. 

Obviously, the resonance wavelength also depends on 
the angle of incidence and it can be pushed down to the 
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FIG. 2: Graphene optical switch based on resonant 
Fabry-Perot transmission, (a) Fabry-Perot resonator in¬ 
corporating a tunable graphene layer inside the cavity flanked 
by two Bragg mirrors (see inset for a sketch of the period and 
labels for geometrical and optical parameters). (b,c) Normal- 
incidence transmittance (b) and reflectance (c) for different 
levels of doping. The cavity is filled with air, but similar per¬ 
formance is achieved with a glass-filled narrower cavity (see 
supplementary Fig. |9|. 


visible regime (Fig. [l]f), although the maximum trans¬ 
mission decreases towards smaller wavelengths due to 
the gradual involvement of interband transitions in the 
graphene. 

In a realistic device, the two graphene layers of Fig. 
[TJd could be biased with a relative potential difference 
V, so that they will reach a Fermi energy \Ep\ = 
hv f y/V cbn /4g?bn 5 where vp « 10 6 m/s is the Fermi ve¬ 
locity in the carbon layer, while cbn and dp n are the 
static permittivity and thickness of the BN layer. For 
dp n ~ 45 nm, a value of Ep = 1 eV is obtained with po¬ 
tentials ^ 4 V. For the single graphene-layer structures 
noted above (see supplementary Figs. [7]and[8|, doping 
could be introduced through the addition of a transpar¬ 
ent electrode. 


B. Graphene optical switch based upon a 
Fabry-Perot resonator 


The concept of the tunneling structure in Fig. [l] can 
be extrapolated to other types of resonators in which the 
incident field also undergoes a large enhancement at a 
position decorated with graphene. A particularly conve¬ 
nient implementation of this idea is presented in Fig. [2j 
as it allows operating under normal incidence conditions. 
More precisely, we replace the tunneling structure by a 
Fabry-Perot (FP) frequency-selective filter, consisting of 
a cavity flanked by two non-absorbing, nearly perfectly 
reflecting mirrors. In practical devices, one generally uses 
Bragg mirrors such as those sketched in Fig. (2^i, which 
are easy to fabricate by multilayer deposition. We con¬ 
sider a separation between the FP mirrors that produces 
a single resonant transmission peak in the 730 — 748 nm 
spectral region. At resonance, light is trapped inside the 
cavity, so it makes many passes through it before escap¬ 
ing, thus generating a large field enhancement at several 
interference nodes. We place the graphene at one of those 
nodes. A similar strategy has been recently followed for 
all-optical nonlinear NIR light modulation [38]. An in¬ 
terplay between absorption (imaginary part of the sus¬ 
ceptibility) and polarization (real part) in the graphene 
is then taking place, leading to large (but not totally 
complementary) modulations in reflection and transmis¬ 
sion, similar to those discussed above for the tunneling 
device. Incidentally, similar performance is obtained by 
filling the cavity with glass and reducing its size (see 
supplementary Fig. [9|, thus configuring a more robust 
structure. We have also verified (results not shown) that 
further reduction of the cavity leads to a ID photonic 
crystal that exhibits a normal-incidence gap, in which a 
localized optical mode exists due to the addition of an 
impurity (i.e., the cavity itself); the graphene can then 
couple to this localized mode to produce an even more 
compact light modulator. We find it interesting that the 
cavity is unaffected if the graphene is placed at an antin¬ 
ode of the interference standing wave inside the cavity 
(see supplementary Fig. 10), as this opens the possibility 
of using an optically inactive graphene layer located at 
an antinode and serving as a gate with which to dope the 
other graphene layer placed at a node. 


C. Enhanced graphene optical absorption and 
switching by coupling to Mie cavities 


Figure [3^i represents the change in the absorption cross 
section undergone by a layer of undoped graphene when 
we place a silicon sphere (e = 12) in its vicinity. These 
types of silicon colloids have been recently synthesized 
and used as excellent photonic cavities m • The increase 
in absorption cross section £cr abs remains a small fraction 
of the extinction produced by the sphere in this configu¬ 
ration (e.g., 6.1% and 2.7% for the Mie modes labeled A 
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FIG. 3: Graphene absorption enhancement by cou¬ 
pling to Mie resonances, (a) Absorption cross section nor¬ 
malized to the projected sphere area ( ttR 2 ), estimated for the 
silicon-sphere/undoped-graphene system shown in the inset 
using Eq. |l]) and Mie theory. We plot the increase in absorp¬ 
tion due to the presence of the sphere. The silicon/graphene 
separation is d — R/ 150. The upper scale corresponds to a 
sphere radius R — 300 nm and d = 2nm. The incident elec¬ 
tric field is along the x direction, (b) Parallel electric-field 
intensity enhancement (\E X \ 2 + \E y \ 2 )/\E 0 \ 2 at the graphene 
plane for the two Mie resonances labeled A and B in (a). The 
quality factors Q of these resonances are also indicated in (a). 


and B in Fig. [ 3 J 1 ), so we approximate it as 

<5(7 abs « 7T a J dxdy |E||/£;o| 2 , (1) 

where E\\ is the parallel component of the electric field 
scattered by the sphere alone, Eq is the incident field, 
and we integrate over the graphene plane. The field 
E ll is obtained from Mie theory m- This approximate 
method yields similar results as the change in elastic 
(dark-field) scattering due to doping, calculated from a 
rigorous modal expansion for the sphere-graphene system 
(see supplementary Fig. 11). In Fig. [ 3 ] the cross section 
is normalized to the projected area of the sphere nR 2 
and the wavelength is normalized to the sphere radius 
R, so that this plot is independent of R, apart from the 
relatively small variations of the permittivity of silicon 
over the NIR. Despite the subwavelength size of the par¬ 
ticle, its high e allows it to trap light within Mie modes 
of high quality factor (Q « 193 and 49 in modes A and 
B, see Fig. If), giving rise to large local enhancements 
of the near-field intensity at the plane of the graphene 
(see Fig. [3 ]d). This in turn boots the absorption, which 
takes remarkably large values, with a peak increase in 
cross section reaching ^ 40% of the projected area of the 
sphere. Interestingly, the spatial distribution of absorp¬ 
tion (proportional to the intensity plotted in Fig. ^jp) is 
strongly confined to the near-contact region, which could 
be exploited for engineering the spatial distribution of op¬ 
tically induced heat deposition, as well as for controlling 
the graphene electron-gas ultrafast dynamics before re¬ 
laxation and thermalization of the absorbed energy takes 
place. 

Because the maximum value of £cr abs produced by a 
single silicon sphere is comparable to its projected area, 
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FIG. 4: Tunable absorption in graphene decorated 
with a 2D array of Mie resonators, (a) Geometry 
and parameters of a triangular array of silicon spheres near 
graphene, (b) Normal-incidence transmission through the 
sphere array without graphene for different lattice periods 
P. The wavelength is shown both normalized to the sphere 
radius R (lower scale) and for R — 300nm (upper scale), 
(c) Absorbance of the array when it is placed near undoped 
graphene (silicon-carbon distance d = R/150) under normal 
incidence. The lattice period is P — 800 nm. We compare two 
different approximations for the graphene conductivity (full 
random-phase approximation [4TJ [42! (RPA) and local-RPA 
nsi) with a semi-analytical model, as discussed in the main 
text. 


we expect to obtain unity-order changes in the absorp¬ 
tion when the graphene is decorated by a periodic array. 
This is illustrated in Fig. [4j where we concentrate on the 
spectral region around the rightmost Mie mode of Fig. 
If (labeled B). We consider the silicon spheres to be ar¬ 
ranged in a triangular lattice, which we simulate using 
a layer-KKR approach [43| (see Methods). Interestingly, 
there is strong interaction between the particles for the 
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lattice spacings P under consideration, which can be in¬ 
tuitively quantified from the fact that the extinction cross 
section of the sphere equals the area of a circle of diameter 
1.75 /im. The transmission of the particle array experi¬ 
ences dramatic spectral variations as P is changed, even¬ 
tually generating a narrow transmission peak, which is 
relatively close, but not on top of the lowest-order Wood 
anomaly, occurring when the wavelength is equal to the 
period at normal incidence; we thus attribute this feature 
to the interaction between Mie modes of the spheres, as 
the wavelength is close (but not right on) a lattice res¬ 
onance that narrows the resulting spectral feature (see 
supplementary Fig. 12). A similar mechanism leading to 
sharp, narrow asymmetric resonances has already been 
described in the context of cavity-waveguide coupling 
M- The absorbance associated with this narrow peak is 
boosted, approaching 50% with undoped graphene (Fig. 
Eh whereas doped graphene shows comparatively negli¬ 
gible absorbance (not shown). 


In the layer-KKR simulation method [43 , the homo¬ 
geneous graphene film enters through its reflection and 
transmission coefficients for different diffracted orders 
(i.e., a collection of propagating and evanescent plane 
waves, each of them corresponding to a fixed value of 
the parallel wave vector). This allows us to use the 
full random-phase approximation (RPA) conductivity 
gu 021 cr(&||, uj), which includes nonlocal effects asso¬ 
ciated with finite parallel wave vectors fey correspond¬ 
ing to those diffracted orders. Because the size of the 
spheres and the lattice periods under consideration are 
large compared with vp/cu (i.e., the ratio of the graphene 
Fermi velocity to the light frequency, ^0.8 nm for a wave¬ 
length of 1.5 /im), nonlocal effects are in fact negligible, 
which explains why we obtain the same results on the 
scale of the figure by just using the value a = e 2 /eh 
for the conductivity in undoped graphene instead of the 
full RPA. The same argument explains why plasmons are 
not excited here in doped graphene. Additionally, we ob¬ 
tain very similar results with the semi-analytical model 
of Eq. 0 (except very close to the resonance), which 
we apply by averaging the parallel electric-field intensity 
enhancement over a unit cell. The intensity in the semi- 
analytical model is calculated without the graphene, just 
to provide insight into the absorption process. However, 
when we calculate it including the carbon layer, the ab¬ 
sorbance A predicted by Eq. 0 cannot be told apart on 
the scale of the figure from the one given by the far-held 
transmittance and reflectance (i.e., A = 1 — T — R), thus 
corroborating the numerical accuracy of our calculations. 


For tutorial purposes, the above discussion on the 
coupling to Mie resonances is based on self-standing 
graphene, but qualitatively similar conclusions are ob¬ 
tained when examining graphene supported on a glass 
substrate (see supplementary Fig. 12). 


D. Enhanced graphene optical absorption and 
switching by coupling to lattice resonances 


We now discuss the absorption enhancement produced 
by lattice resonances, for which strong scatterers such 
as metallic particles are preferable. Although metals in¬ 
troduce additional losses, their absorbance is relatively 
small in the NIR, so graphene can still make a big dif¬ 
ference. This is corroborated in Fig. [5j where we con¬ 
sider a graphene layer decorated with a 2D square ar¬ 
ray of gold spheres surrounded by silica for different val¬ 
ues of the lattice spacing P. The transmission (Fig. J5j)) 
and reflection (Fig. spectra of these structures ex¬ 
hibit sharp features emerging near the Wood anomaly 
condition (i.e., when the wavelength in the surrounding 
dielectric is close to the period, or equivalently, when a 
diffraction order becomes grazing), which can be easily 
understood in terms of lattice resonances [45J 0B]. As 
the period is increased, these features move to the red, 
where the metal is less lossy, and consequentally, the reso¬ 
nances become narrower. The additional absorption pro¬ 
duced by the undoped graphene then becomes more no¬ 
ticeable (see supplementary Fig. 13), eventually causing 
a decrease in peak transmittance of ^ 60%, accompanied 
by a 28-fold reduction in reflectance. 


III. CONCLUSIONS AND OUTLOOK 

In conclusion, monolayer graphene can be used to 
produce unity-order changes in the transmission, reflec¬ 
tion, and absorption of light down to the vis-NIR do¬ 
main when comparing its electrically doped and undoped 
states, considering realistic levels of doping (Ep ~ 1 eV) 
that are currently attainable through electrical gating. 
It should be emphasized that the calculations here pre¬ 
sented for geometries containing only graphene and di¬ 
electrics are scalable, so that the main requirement is 
that the graphene can be switched back and forth be¬ 
tween Ep = 0 and Ep > E p / 2, where E v is the photon 
energy under operation; provided this condition is sat¬ 
isfied, all geometrical lengths and the light wavelength 
can be scaled by a common factor, leading to the same 
values for the transmission and absorption. For example, 
the modulation at 700 nm wavelength predicted in Fig. [l] 
with doping up to Ep = 1.1 eV can be also extrapolated 
to the 1550 nm telecom wavelength with doping up to 
Ep = 0.5 eV after scaling all lengths by a factor of ^ 2. 

Interestingly, we find that undecorated graphene in a 
planar multilayer dielectric structure can modulate trans¬ 
mission near the point of resonant tunneling under total 
internal reflection, with absolute changes exceeding 90% 
and an extinction ratio > 15 dB. Similar levels of mod¬ 
ulation are found for graphene placed inside a realistic 
Fabry-Perot cavity. Large vis-NIR modulation depths 
are also predicted in a graphene layer decorated with pe¬ 
riodic arrays of silicon or gold particles. Obviously, the 
depth of modulation is reduced by coupling to impuri- 
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FIG. 5: Enhanced tunable graphene absorption by coupling to lattice resonances in 2D metal particle arrays. 

(a) Square array of gold spheres (radius R) placed above graphene (2nm gold-to-carbon separation). The entire system is 
assumed to be embedded in silica (e = 2.25). (b,c) Normal-incidence transmission (b) and reflection (c) spectra for R — 80 nm 
and different lattice periods P with either doped (broken curves, Ef = 1 eV) or undoped (solid curves) graphene. The spectra 
are dominated by lattice resonances occurring near a free-space light wavelength A ~ PyT. (d) Peak wavelength with doped 
graphene (right scale) and transmission at that wavelength with either doped or undoped graphene (left scale) as a function of 
gold sphere radius for a period P = 500 nm. (e) Same as (d) for silver particles. 


ties in low-quality graphene, where optical losses can be 
still significant under high doping. Nonetheless, we find a 
substantial degree of modulation even in the presence of 
large residual absorption (e.g., ^ 50% modulation in the 
device of Fig. [l] when the residual optical loss amounts to 
14% of the ideal absorption-free highly doped material, 
as shown in Fig. [l^). 

The mechanisms here considered for light modulation 
by graphene can be integrated in devices spanning only 
a few square microns in size, so they require a relatively 
small amount of doping charge to operate. We thus an¬ 
ticipate that these systems will be able to modulate vis- 
NIR light at high speeds with a minute consumption of 
power, typical of capacitive devices. This is an advantage 
with respect to alternative commutation devices based on 
quantum-wells m and phase-change materials [58] . For 
example, we envision an integrated commutation device 
operating over an area A = 50 x 50 fim 2 (i.e., covering 
a customary optical beam size), for which we estimate 
a capacitance C = Ae/Aird ~ 0.3 pF, where we consider 
e = 4 (DC silica) and a gate separation d = 300 nm (no¬ 


tice that there is great flexibility in the choice of d in some 
of our devices). The time response is then limited by 
the sheet resistance of the graphene layer (^ 100 !!/□), 
giving an overall cutoff frequency for the electrical band¬ 
width of 1/2ttRC ~ 5 GHz, while the optical limit for 
the electrical modulation of the photonic response (i.e., 
the effect related to the decay time of the resonance) 
renders a larger cutoff (c/2LQ ~ 150 GHz for a cavity 
length L ~ 1 fim and a quality factor Q ~ 10 3 ). The 
large electro-optical response of graphene combined with 
its small volume are thus ideal attributes for the design 
of fast optical modulators and switches operating in the 
vis-NIR, which can benefit from the coupling to optical 
resonators such as those explored in the present work. 
In particular, the planar structures presented in Figs, [l] 
and[2j which rely on unstructured graphene, provide rel¬ 
atively affordable designs that are appealing for micro 
integration and mass production. 






















METHODS 

Schrodinger equation for TE polarized light in 
a planar multilayer structure. For our purpose, it 
is convenient to write the wave equation for the electric 
field EasVxVxE — & 2 eE = 0, where k = uj/c is the 
free-space light wave vector and e is the frequency and 
position dependent local dielectric function. We consider 
a structure formed by several planar layers (see Fig. [TJ, 
illuminated under TE polarization (i.e., with the elec¬ 
tric field parallel to the planes), and oriented with the 
z axis perpendicular to the interfaces, so that the spa¬ 
tial dependence of the dielectric function is only through 
z (i.e., e(z)). Then, we can write the electric field as 
E(r) = 'ip(z)e lk w y it, where fey = ky/e[ sin 0{ is the paral¬ 
lel wave vector component, which is determined by the 
angle of incidence 6{ at the near medium of permittivity 
€[. With this form of the electric field, the wave equation 
reduces to 

+ v { Z )^ z )=mz\ 

where we have defined the equivalent potential V(z) = 
[1 — e(z)]k 2 /2 and the effective energy E = (k 2 — k 2 )/ 2. 
Interestingly, metals (e < 0) and dielectrics (e > 1) pro¬ 
duce repulsive (V > 0) and attractive (V < 0) potentials, 
respectively, in this equivalent Schrodinger model. 

Graphene conductivity. We model graphene 
through its 2D conductivity. For doped graphene, we 
use a previously reported local-RPA model 03 023] , con¬ 
veniently corrected to account for finite temperature T 
in both intra- and interband transitions m- More pre¬ 
cisely, 

e 2 i 

^ ^ 7T H 2 (i0 + iT -1 ) 

dE _ /# ~ f-E _ 

1 — AE 2 /[h 2 (uj + ir -1 ) 2 ]J 5 

where Ep is the Fermi energy, 

Ep = Ep + 2 k B T log (l + e - BF / feBT ) 

effectively accounts for thermal corrections in the doping 
level, and /# = 1/[1 + e ( E ~ E F)/kBT j is the Fermi-Dirac 
distribution function. For undoped graphene, the above 
expression converges to the well-known limit a — e 2 /Eh. 
We further account for finite wave-vector effects (non- 
locality) through the full-RPA model [HJ 02] , in which 
<j depends on fey and cj, but we find those effects to be 
negligible (see Fig. [4|, as expected from the large mis¬ 
match between the Landau damping range uf/cj at vis- 
NIR frequencies uo and both the light wavelength and the 
distances involved in the structures under consideration. 
Throughout this work, we take T = 300 K and assume 
an inelastic decay time given by the Drude model EDE21 
(i.e., r = fiE-p/eVp) with an impurity limited DC mobil¬ 
ity /i = 2000 cm 2 / (V s). 


Multilayer structure simulation. The transmis¬ 
sion and field enhancement of planar multilayers are ob¬ 
tained through a standard transfer matrix approach. In 
particular, we use the reflection and transmission coef¬ 
ficients for a plane wave of parallel wave vector fcy in¬ 
cident from medium 1 on a graphene layer of conduc¬ 
tivity a placed at the interface with another medium 2, 
which upon direct solution of Maxwell’s equations for 
the s (TE) polarization under consideration are found 
to be m r i 2 ~ (kzi ~ k Z 2 + 9s) /{k z i + k z 2 + g s ) and 
t s 12 = 2k z i/(k z i + k z 2 E g 8 ), respectiv ely, where g s = 

Atujuj/c 2 and k z j = yjk 2 ej — k 2 + i0+. For complete¬ 
ness, we give the coefficients for p (TM) polarization: 
r i 2 = ( e 2 ^i - e x k z2 + 9 P )/ (e 2 k zl + e x k z2 + g p ) and t p 12 = 
Zy/ei&^kzi/(e 2 k z i-\-eik z2 +g p ), where g p = A7rak zl k z2 /uj. 
Incidentally, the sign of the square root is chosen to yield 
positive real values. These expressions also describe the 
coefficients of interfaces without graphene, simply by tak¬ 
ing g s = 9 P = 0. 

For tunneling transmission (Fig. [I]), the resonant wave¬ 
length of maximum transmission A res is only slightly 
changed from the central waveguide Fabry-Perot reso¬ 
nance condition, k z2 d+tp = ALr, where d is the waveguide 
thickness, N (= 0 under the conditions of this work) is 
the order of the resonance, and p = arg{r 21 } (we choose 
media 1 and 2 right outside and inside the waveguide, 
respectively). To linear order in cr, we find 

-2Trk z2 d ( 87r(fc/c)cos(v?/2)Im{<r}\ 

TEf+TSF )■ (2) 

Incidentally, a has units of velocity in CGS, so this ex¬ 
pression is dimensionally correct. 

Particle arrays. We use the layer-KKR method to 
simulate periodic particle arrays near planar interfaces 
m- This method relies on an expansion of the electro¬ 
magnetic field in terms of spherical vector waves around 
the particles and plane waves near the graphene. The 
scattering by the spheres then involves multiplication by 
Mie coefficients, whereas the graphene enters through its 
reflection coefficients (see above). Plane and spherical 
waves are analytically transformed into each other, giv¬ 
ing rise to a self-consistent system of equations projected 
on the coefficients of the sphere multipoles. Translational 
lattice symmetry is used to reduce the number of plane 
waves to those of a discrete set corresponding to different 
diffraction orders (i.e., two waves of orthogonal polariza¬ 
tions for each reciprocal lattice vector). We achieve con¬ 
vergence with ~ 100 such waves and neglecting sphere 
multipoles of orbital angular momentum number above 
7. This method directly yields the reflection, transmis¬ 
sion, and absorption coefficients used to produce Figs. [4] 
and [5] for periodic particle arrays near planar interfaces 
including graphene. 
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FIG. 6: (Supplementary figure.) Reflectance (a), reflection phase (b), and transmission phase (c) under the same conditions as 
in Fig. le of the main paper. 
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FIG. 7: (Supplementary figure .) Transmittance (b) and reflectance (c) for a structure similar to that of Fig. 1 of the main 
paper, but containing only one graphene layer (see (a)), under the same conditions of light incidence. 
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FIG. 8: (Supplementary figure.) Reflectance (b) and absorbance (c) of a similar structure as in Fig. 1 of the main paper, but 
containing only one graphene layer and without any out-coupling BF11 medium on the right-hand side of the structure. The 
conditions of light incidence are the same as in Fig. 1. 
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FIG. 9: (Supplementary figure .) Performance of a Fabry-Perot cavity similar to that of Fig. 2 of the main paper, but filled 
with glass and designed to operate in the same spectral region using modified geometrical parameters. 
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FIG. 10: (Supplementary figure.) Electric field intensity enhancement relative to the incident intensity inside the Fabry-Perot 
cavity considered in Fig. 2 of the main paper, calculated at the 738 nm resonance wavelength in the absence of graphene. 
The addition of a second graphene layer at an antinode (rightmost graphene layer in this plot) produces exactly the same 
transmission and reflection spectra as in Fig. 2, regardless the doping state of the extra layer. The width of the cavity is 800 nm 
and other geometrical parameters are the same as in Fig. 2a. 
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FIG. 11: (Supplementary figure.) Graphene-doping-induced change is the dark-field scattering cross section (i.e., integrated 
over scattering directions other than the specular reflection or direct forward transmission) of a silicon sphere under the same 
conditions as in Fig. 3a. The curve represents the difference in angle-integrated elastic cross section when the graphene is 
undoped or doped to Ep = 1 eV, calculated using a modal expansion described elsewhere m- 
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FIG. 12: (Supplementary figure.) (a) Normal-incidence (fey = 0) absorption spectra for a triangular lattice of silicon spheres 
(radius R — 300 nm and lattice period P = 800 nm) placed on top of a graphene sheet (silicon-carbon separation distance 
d— 2 nm) when the graphene is supported on a silica substrate, (b) Same as (a) for a square lattice, (c) Dispersion diagram 
of the triangular silicon-sphere lattice without graphene in the Mie resonance region under consideration. The white vertical 
segment in (c) indicates the spectral range in (a), dominated by a sphere Mie mode that is crossed by a lattice resonances at 
finite k *. The lattice resonance produces a narrowing of the Mie mode. 
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FIG. 13: Absorbance spectra under the same conditions as in Fig. 5b,c of the main paper. 














